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Abstract. We present the results of high S/N long-slit spectroscopy with the Multiple Mirror and the SAO 6-m 
telescopes, optical imaging with the Wise 1-m telescope and H l observations with the Nangay Radio Telescope of 
the very metal-deficient (12-|-log(0/H)=7.64) luminous (Mb=-18™1) blue compact galaxy (BCG) HS 0837-1-4717. 
' The blue bump near A4650, characteristic of WR stars is detected in the central supergiant H ll region, as well 

I as the barely seen red bump at A5808. The derived number of WR stars in the region of the current starburst is 

. ~1000. Evidence for fast motions in this region is also seen as broad, low-contrast components in the Ha, H/3 and 

' strong [O III] lines AA4959,5007. While the extinction of the narrow emission lines from the supergiant H ll region 

is low, the very large Balmer decrement of the broad components suggests that the part of current starburst is 
highly obscured by dust. Abundance ratios X/0 for X=Ne, Ar, S, Fe and CI in the supergiant Hll region are in 
O ' good agreement with the mean values of other very metal-deficient BCGs. Nitrogen, however, is overabundant by 

a factor of ~6. This implies an unusually efficient N enrichment in HS 0837-1-4717, and probably, a non-typical 
C/3 ' evolution scenario. The Ha-line Position- Velocity (P-V) diagrams for directions approximately along the major 

^ ' and minor axes reveal disturbed motions of the ionized gas, mainly in peripheral regions. The SW part of the major 

axis P-V diagram looks like a rotation curve, with the velocity amplitude Kot ~50-70 km s~^ at r=4.3 kpc. Its 
1 J - NE part displays, however, strong deviations, indicating either counter-rotation, or a strong outflow/supershell. If 

r> I it is considered as indicating a shell-like feature its velocity amplitude of ~70 km s~^ (relative to the extrapolated 

^ . rotation curve), and the apparent extent of ~4" (3.3 kpc) imply a dynamical age of ~14 Myr and the full energetic 

equivalent of ~2.3xl0* SNe. The latter indicates continuing starbursts during at least this time interval. The 
long-slit spectra reveal a complex morphology for this galaxy. It consists of two compact regions at a distance of 
~2 kpc. Their continuum flux differs by a factor of four. The brightest one is related to the current starburst with 
the age of ~3.7 Myr. The slightly redder fainter component could be an older starburst (~25 Myr). The Wise 1-m 
telescope UBVR integrated photometry reveals a high optical luminosity for this BCG, and the unusual (B — V^) 
and {V ~R) colours. The morphology of HS 0837-)-4717 is highly disturbed, with two small tails emerging to NNW 
and SSE. Such a disturbed overall morphology, a " double- nucleus" structure, significantly disturbed velocities of 
ionized gas, together with the very high power of the starburst suggests a possible explanation of the object as 
a recent merger. We compare the properties of this BCG and of similar objects known in the literature, and 
conclude that their high nitrogen excess is most probably related to the short phase of a powerful starburst when 
many WR stars contribute to the enrichment of ISM. 

Key words, galaxies: starburst - galaxies: abundances - galaxies: interactions - galaxies: evolution - stars: Wolf- 
Rayet - galaxies: individual (HS 0837-1-4717) 



1. Introduction 



Most low-mass gas-rich galaxies have low metallicities. 
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of Kunth & Ostlin l7l)Ul)|l . Only for a very small fraction 
of galaxies in the Local Universe is the metallicity ex- 
tremely low, in the range of Zq/50 to Zq/20. This range is 
more characteristic of high-redshift damped Lya absorp- 
tion systems which could be young galaxies. The study of 
the extremely metal-deficient (XMD) galaxies allows one 
to test many theoretical ideas on massive star formation 
and evolution, and models of galaxy evolution in a very 
low metallicity environment. 

In particular, the abundance patterns of a-elements 
and iron, that is, the ratios X/0 (with X=Ne, S, Ar, Fe), 
allow one to check models of stellar nucleosynthesis. For 
metal-poor galaxies they were found to be remarkably con- 
stant over a wide range of 0/H (e.g., Izotov & Thuan [T999l 
hereafter IT99). This implies that all these elements are 
primary and are produced in the same massive stars. The 
ratio N/0 was also found to be fairly constant for XMD 
BCGs (IT99), implying mainly primary production of ni- 
trogen in very low metallicity environment. 

Wolf-Rayet (WR) stars - a specific very short phase 
of massive star evolution, characterized by strong mass 
outflow, and enriched mainly in N and C, are often seen 
in high S/N spectra of starburst galaxies (e.g., Kunth & 
Joubert fnSS ! Izotov et ak HMHI Guseva et al HUU^ . These 
stars are detected through the broad emission features of 
the so called "blue bump" near A4650 A, and more rarely 
in the red, near A5808 A. Models of stellar evolution pre- 
dict that the number of WR stars in a young star cluster 
is a sensitive function of metallicity, IMF and starburst 
age (e.g., Schaerer & Vacca lTMHl hereafter SV98). 

WR star outflows could produce a significant N over- 
abundance at the locations of young starbursts. However, 
the analysis of a large sample of Hii galaxies by 
Kobulnicky & Skillman t|1996) showed that galaxies, with 
and without strong WR features in their integrated spec- 
tra, show identical N/0 ratios. This is probably explained 
by nitrogen outflow in the hot phase. Therefore, only sig- 
nificantly after its ejection could the nitrogen be cooled 
and mixed into the optically observed ionized gas. 

Significant nitrogen overabundance is detected in a few 
starburst galaxies. The most prominent one is Mkn 996 
(Thuan et al. HlMjl . with 12-)-log(O/H)=8.0. It is a non- 
typical BCG due to the large number of WR stars and to 
a powerful outflow from the compact central SF burst. Its 
N/0 in the region with diameter of 0.6 kpc is 4 to 25 times 
higher than the typical value for low-metallicity BCGs of 
^1/40. The morphology of Mkn 996 implies that likely it 
is a remnant of a recent merger. The authors argue that 
the observed large N/0 ratio is directly related to the 
powerful outflows of numerous WR stars. 

Another well known case is a factor of three nitrogen 
overabundance, found in two compact regions (d ~20 pc) 
near the central starburst in the nearby (D ~4 Mpc) dwarf 
WR galaxy NGC 5253 (e.g., Kobulnicky et al. lTMTl and 
the summary of previous results therein). 

Only three galaxies: UM 420, Mkn 1089, and UM 448, 
all luminous and with detected WR lines (Guseva et al. 
IM)0,1 . from the sample of 50 BCGs in IT99, show nitrogen 



overabundance of a factor of «3 relative to the mean N/0 
value of 1/30 for their range of O/H). Finally, the galaxy 
Haro 11 from the sample of luminous BCGs (Bergvall & 
Ostlin IT()()2|I . also with strong WR bump, shows an N/0 
ratio that is 6 times higher than for the main BCG group. 

In this paper we present the results of a comprehen- 
sive study of HS 0837-^4717, found in the Hamburg/SAO 
survey (Pustilnik et al. I1999|l . The very low metallicity 
of HS 0837-^4717 {Z kZq/20) was first derived from the 
observations with the SAO RAS 6-m telescope (BTA) in 
1996 (Kniazev et al. SOOOa). Classified as an XMD BCG, 
it is unusually luminous for its metallicity. This galaxy 
shows WR bumps and broad components of strong emis- 
sion lines of hydrogen and oxygen. This is the next most 
nitrogen overabundant galaxy known after Mkn 996. But, 
in contrast to the latter, its optical spectrum does not look 
atypical for BCGs, and its very large nitrogen overabun- 
dance is seen on an otherwise more or less typical BCG 
background. 

In Sect. we describe the observations and their re- 
duction. In Sect. 13 all the results from the reduction and 
preliminary analysis are presented, including the heavy 
element abundances, the WR and broad components of 
strong lines, the ionized gas kinematics and the morphol- 
ogy and photometry. In Sect. ^ we discuss the obtained 
results and derived parameters of this galaxy, compare its 
properties with those of other luminous and XMD BCGs, 
and draw conclusions. The adopted distance to the galaxy 
is 170.5 Mpc; respective scale is 827 pc per arcsecond. 

2. Observations and data reduction 

2.1. Spectrophotometry and kinematic data 

Spectrophotometric observations of HS 08374-4717 with 
high S/N ratio were obtained with the MMT on May 
1, 1997. The observations were performed with the Blue 
Channel of the MMT Spectrograph using a highly- 
optimized Loral 3072x1024 CCD detector. A l'.'5xl80" 
slit was used along with a 500 grooves nim^^ grating used 
in the first order and with the L-38 second order block- 
ing filter. This yields a spatial scale along the slit of 0'.'3 
pixel^^, a scale perpendicular to the slit of 1.9 A pixel"^, a 
spectral range 3600-8300 A and a spectral resolution of ~ 
7 A (FWHM). For these observations the CCD rows were 
binned by a factor of 2, yielding a final spatial sampling 
of 0'.'6 pixel^^. The observations cover the full spectral 
range in a single frame which contains all the lines of in- 
terest. Furthermore, the spectra have sufficient spectral 
resolution to separate [O ill] A4363 from the nearby H7 
and to distinguish between the narrow nebular and the 
broad WR emission lines. The total exposure time was 40 
minutes, split into 2 equal subexposures, to allow for more 
effective cosmic ray removal. Both exposures were taken at 
an airmass of ~1.13. The slit was oriented in the direction 
with the position angle PA = 16°, along the major axis 
of the brightest region of the galaxy (Fig. ^ . The spec- 
trophotometric standard star HZ 44 (Massey et al. I1988|l 
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Telescope 


Date 


Exposure 


Wavelength 


Dispersion 


Seeing 


Airmass 


PA 






time [s] 


Range [A] 


[A/pixel] 


[arcsec] 




[degree] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


MMT 


01.05.1997 


2x1200 


3600 - 8300 


1.9 


1.6 


1.13 


16 


BTA 


19.02.2002 


2x1800 


3700 - 6100 


2.4 


2.2 


1.27 


44 


BTA 


19.02.2002 


2x1800 


6000 - 7200 


1.2 


1.3 


1.15 


44,130 



T 1 1 1 1 1 1 I 1 1 i^y I r 




-5 5 

Position (arcsec) 

Fig. 1. Morphology of HS 0837+4717 as seen in the Wise 
V^-band image. North is up, East is to the left. The po- 
sitions of MMT (PA=16°) and BTA (PA=44° and 130°) 
long slits are shown. The isophotes sample the image at 
the brightness levels of 11, 20, 30, 50, 75, 100, 150, 200 
and 220 units. The minimal level corresponds to ^1.5(7 of 
the noise. See Sect. 13.51 

was observed for flux calibration. Reference spectra of He- 
Ne-Ar were obtained before and after each observation to 
provide wavelength calibration. 

Observations with the BTA were conducted on 
February 19, 2002 with the Long-slit spectrograph 
(Afanasiev et al. I1995|l at the prime focus of telescope. 
A CCD detector Photometries PM1024 with 24x24 ^m 
pixel was used for data acquisition. For abundance study 
a grating with 651 grooves mm~^ was used in the first 
order, giving 2.38 A pixel-i (FWHM=7 A) and allow- 
ing the simultaneous coverage of the range 3700-6100 A. 
Two 30-min exposure spectra were averaged to yield the 
final spectrum. PA=44° was selected to be close to the 
apparent major axis of the object's fainter isophotes. To 
complement the study of emission lines in the blue, and to 
study the ionized gas kinematics, we obtained two spectra 
in the red with a grating with 1302 grooves mm~^, giving 
1.2 A pixel-\ and FWHM=3.8 A. These spectra covered 



the range 6000-7200 A, and were obtained also with the 
slit direction along PA=44°, and in the nearly perpendic- 
ular direction of PA=130°. The long sht of 2"xl20" was 
used for both observations. 

The reduction of the spectral observations was carried 
out at SAO using the IRAF^ and the MIDAS^ software 
packages. It included bias subtraction, cosmic ray removal, 
and flat-field correction using exposures of a quartz in- 
candescent lamp. After wavelength mapping, night sky 
background subtraction, and correcting for atmospheric 
extinction, each frame was calibrated to absolute fluxes. 
One-dimensional spectra were extracted by summing up, 
without weighting, different numbers of rows along the 
slit depending on exact region of interest. The MMT spec- 
trum of the brightest region of HS 0837+4717 in the range 
3600 A < A < 7300 A for the aperture l'.'5x3'.'0 and the 
BTA spectra in the blue and red, extracted with the aper- 
ture 2"x2'.'8, are shown in Fig. [3 

2.2. Photometry data 

Photometric observation of HS 0837+4717 were performed 
with the Wise observatory 1-m telescope on nights 27 and 
28.05.1999. in U,B,V and Rc bands. The images were 
exposed through the standard filters onto a Tektronics 
1024x1024 pixels CCD and Ilf9xllf9 field-of-view (at 
0.696 arcsec pixel"^). During the first night, the obser- 
vations were performed at airmasses of 1.69 till 2.3, be- 
ginning from the [/-band, and ending with the i?c-band. 
The exposure time was 30 min for U, and 15 min for the 
others filters. For the second night, the observations were 
collected at the airmasses from 1.54 to 2.20, with the same 
exposure times. Photometric standards from the fields of 
Landolt (1992J were observed at an airmass range close to 
that of the studied galaxy. The seeing during both nights 
was ~(2-2.5)". 

The reduction of photometry was performed in IRAF, 
following the standard pipe-line, including the de-biasing, 
cosmic ray removal and flat-fielding. Aperture photome- 

^ IRAF: the Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for 
Research in Astronomy, In. (AURA) under cooperative agree- 
ment with the National Science Foundation (NSF). 

^ MIDAS is an acronym for the European Southern 
Observatory package - Munich Image Data Analysis System. 
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Fig. 2. Upper panel: The MMT 2D spectrum of HS 0837+4717 with PA=16°, plotted with the contrast level adjusted 
so as to show the distribution of the continuum intensity along the slit. SSW is up. While the intensity of strong 
lines decreases symmetrically relative to the continuum brightness peak, the continuum intensity at the lower level is 
skewed to the SSW direction, indicating a slightly non-central position of the starburst. Middle panel: De-redshifted 
MMT ID spectrum of the bright Hii region in HS 0837+4717, extracted with the aperture of r.'5x3'.'0. Bottom panel: 
De-redshifted ID blue and red BTA spectra. Note the good consistency between the MMT and BTA spectra. 
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Ao(A) Ion 


MMT (l'.'5x3'.'0) 






BTA(2" 


x2'.'8) 




F(A)/F(H/3) 


7(A)//(H/3) 


F(A)/F(H/3) 


/(A)/7(H/3) 


3727 [0 ii] 


0.380±0.019 


0.402±0, 


.022 





.379±0, 


.021 


0.407±0 


.024 


3835 H9 


0.052±0.005 


0.082±0, 


.009 


0, 


.050±0, 


.005 


0.078±0 


.009 


3868 [Nc III] 


0.418±0.021 


0.436±0, 


.024 





.398±0, 


.023 


0.423±0 


.026 


3889 He I + H8 


0.147±0.011 


0.181±0, 


.015 





.145±0, 


.012 


0.178±0 


.016 


3967 [Ne iii] + H7 


0.269±0.013 


0.304±0. 


.016 


0, 


.257±0. 


.014 


0.295±0, 


.017 


4026 He I 


0.018±0.003 


0.019±0, 


.003 





.018±0, 


.003 


0.019±0 


.003 


4069 [S ii] 


0.010±0.003 


O.OlOiO, 


.003 





.012±0, 


.003 


0.013±0 


.003 


4101 US 


0.230±0.011 


0.260±0, 


.013 


0, 


.221±0, 


.012 


0.255±0, 


.015 


4340 H7 


0.466±0.021 


0.490±0. 


.023 


0, 


.449±0. 


.023 


0.480±0, 


.026 


4363 [0 III] 


0.171±0.008 


0.172±0. 


.009 


0, 


.163±0. 


.010 


0.166±0, 


.010 


4471 He I 


0.047±0.004 


0.047±0, 


.004 





.037±0, 


.003 


0.037±0 


.003 


4658 [Fc III] 


0.004±0.002 


0.004±0, 


.002 





.003±0, 


.001 


0.003±0 


.001 


4686 He 11 


0.019±0.002 


0.019±0, 


.002 


0, 


.019±0, 


.002 


0.019±0 


.002 


4701 [Fe III] 


O.OOliO.OOl 


O.OOliO, 


.001 





.004±0. 


.002 


0.004±0 


.002 


4713 [Ar iv] + He i 


0.023±0.002 


0.023±0, 


,002 





.022±0, 


.002 


0.022±0 


.002 


4740 [Ar iv] 


0.015±0.002 


0.015±0, 


.002 





.015±0, 


.002 


0.015±0 


.002 


4861 H/3 


1.000±0.045 


l.OOOiO, 


.047 


1 


.OOOiO, 


.058 


l.OOOiO 


.059 


4907 [Fe iv] 


0.006±0.002 


0.006±0, 


.002 


0, 


.007±0, 


.003 


0.007±0, 


.003 


4922 He I 


0.012±0.002 


O.OlliO, 


.002 





.014±0, 


.003 


0.013±0 


.003 


4959 [O III] 


2.029±0.089 


1.950±0, 


.089 


1, 


.966±0, 


.099 


1.923±0 


.099 


4989 [Fe iii] 


0.005±0.002 


0.005±0, 


.002 


0, 


.007±0, 


.003 


0.006±0 


.003 


5007 [O III] 


6.089±0.268 


5.834±0, 


.266 


5 


.900±0, 


.300 


5.751±0 


.298 


5041 [Fe iv] 


0.007±0.002 


0.006±0, 


.002 


0, 


.OlOiO. 


.003 


O.OlOiO 


.003 


5047 He I 


0.004±0.002 


0.004±0, 


.002 





.003±0, 


.003 


0.003±0 


.003 


5158 [Fe 11] 


0.006±0.001 


0.005±0, 


.001 





.004±0, 


.002 


0.004±0 


.002 


5199 [N i] 


0.006±0.002 


0.005±0, 


.002 


0, 


.005±0, 


.002 


0.005±0, 


.002 


5271 [Fe III] 


0.009±0.003 


0.008±0, 


.003 


0, 


.012±0. 


.003 


O.OlliO, 


.003 


5518 [CI III] 


0.004±0.001 


0.004±0. 


.001 












5538 [CI III] 


0.003±0.002 


0.003±0, 


.002 












5755 [N 11] 


0.009±0.003 


0.009±0, 


.002 












5876 He I 


0.140±0.009 


0.127±0, 


.009 


0, 


.140±0, 


.008 


0.130±0 


.008 


6300 [O i] 


0.024±0.004 


0.022±0. 


.004 


0, 


.021±0. 


.004 


0.019±0, 


.004 


6312 [S III] 


0.011±0.004 


O.OlOiO. 


.004 


0, 


.016±0. 


.004 


0.014±0, 


.003 


6364 [O i] 


0.008±0.002 


0.007±0, 


.002 





.006±0, 


.002 


0.005±0 


.002 


6548 [N 11] 


0.027±0.002 


0.024±0, 


.002 


0, 


.027±0, 


.003 


0.024±0 


.003 


6563 Ha 


3.126±0.139 


2.756±0, 


.137 


3 


.075±0, 


.156 


2.758±0, 


.154 


6584 [N it] 


0.081±0.012 


0.071±0, 


.011 





.080±0, 


.011 


0.072±0 


.010 


6678 He i 


0.035±0.003 


0.031±0, 


.003 





.034±0, 


.003 


0.030±0 


.003 


6717 [S 11] 


0.043±0.004 


0.038±0, 


.003 


0, 


.046±0, 


.004 


0.041±0, 


.004 


6731 [S 11] 


0.037±0.004 


0.032±0, 


.003 


0, 


.040±0. 


.004 


0.035±0, 


.004 


7065 He i 


0.081±0.005 


0.070±0, 


.005 












7136 [Ar III] 


0.040±0.004 


0.034±0. 


.003 












C(H/3) dex 


0.12±0.06 








0.12±0.07 




EW(abs) A 


3.75±0.41 








3.55±0.42 




F(H/3)'' 


258±3 








313±9 




EW(H/3) A 


230± 6 








222± 9 





" in units of 10 ergs s cm 



try was performed in circular apertures, and the asymp- 
totic values of the growth curves were accepted as the 
final integrated magnitudes. The nights were not perfect; 
some cirrus clouds presumably were present. Therefore, 
the maximal differences, obtained for the magnitudes in 
these 2 nights reach 0™3 (in the V-band). We combined 
the results for both nights in the attempt to obtain more 
reliable magnitudes of this galaxy. We also used as an in- 



dependent control, the BVR magnitudes and colours de- 
rived from the convolution of the MMT spectrum with the 
B, V and R passbands. 
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2.3. H I data 

H I line observations were carried out in February 1998 
with the 300m Nangay^ radio telescope (NRT). The NRT 
has a half-power beam width of 3'.7 (EW) x 22' (NS) at 
the declination Dec. = 0°. 

Since HS 0837-1-4717 had a known optical redshift, we 
split the 1024-channel autocorrelator in two halves and 
used a dual-polarization receiver to increase the S/N ra- 
tio. Each correlator segment covered a 6.4 MHz band- 
width, corresponding to a 1350 km s^^ velocity cover- 
age, and was centered at the frequency corresponding to 
the optical redshift. The channel spacing was 2.6 km s~^ 
before smoothing and the effective resolution after aver- 
aging pairs of adjacent channels and Hanning smoothing 
was 10.6 km s~^. The system temperature of the receiver 
was w 40 K in the horizontal and vertical linear polariza- 
tions. For the target galaxy declination the telescope gain 
was 0.88 K/Jy. The observations were made in the stan- 
dard total power (position switching) mode with 1-mimite 
on-source and 1-minute off-source integrations. The data 
were reduced using the NRT standard packages, written 
by the observatory's staff. The Horizontal and Vertical 
linear polarization spectra were calibrated and processed 
independently, and were averaged together. 

3. Results 

3.1. Chemical abundances 

To derive element abundances of species O, Ne, N, S, 
Ar and CI we use the standard method from AUer 
l|1984|l . and follow the procedure described in detail by 
Pagel et al. H1992|l and Izotov et al. H1994I hereafter 
ITL94). The electron temperature Te is known to be dif- 
ferent in high- and low-ionization Hll regions (Stasihska 
[TM^ . In the following we have chosen to determine 
Te(Olll) from the [O iii]A4363/(A4959+A5G07) ratio us- 
ing the five-level atom model of Oiil, and No(Sii) from 
the [Sli]A6717/A6731 ratio. We then adopt Te(Oiii) for 
derivation of 0^+, Ne^"*" and Ar^+. To derive the elec- 
tron temperature for 0+ ion, we used the relation between 
Te(0 II) and To(0 ill) fitted by ITL94 for the photoionized 
Hii models of Stasihska (1990). re(Oii) has been used to 
derive the 0+, N+ and Fe"*" ionic abundances. For Ar^"*" 
and S^"*" we have used an electronic temperature inter- 
mediate between Te(Oli) and Te(Oiii) following the pre- 
scription of Garnett ()1992|l . 

Total element abundances have been derived after cor- 
rection for unobserved stages of ionization as described by 
ITL94. For the case apphcable to HS 0837-1-4717, when the 
strong nebular Hen A4686 emission line is present, this 
implies that a non-negligible amount of O'^^ ion should 

^ The Nangay Radioastronomy Station is part of the Paris 
Observatory and is operated by the Ministere de I'Education 
Nationale and Institut des Sciences de I'Univers of the Centre 
National de la Recherche Scientifique. 
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Fig. 3. Gaussian multicomponent fitting (solid lines) of 
the BTA spectra of HS 0837-^4717 (dotted lines), top 
panel: region of the blue bump; middle panel: region of 
H/3 and [Oiii] AA4959,5007; bottom panel: region of Ha, 
[Nil] AA6548,6583. 

present. Its abundance is derived from the relation: 

0^+ = He2+/He+(0+ + 0^+). (1) 
Then, the total oxygen abundance is equal to 

= 0+ + 0^+ + 0^+. (2) 
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Fig. 4. Distribution of log(N/0) ratio versus 
12-|-log(0/H) for the BCG sample from Izotov & 
Thuan 0^291 with the value for HS 0837+4717 shown. 
Data for Haro 11 and Mkn 996 are from Bergvall & Ostlin 
Pin^ and Thuan et al. l(TMS|l . respectively. The three 
galaxies with log(N/0) in the range -1.1 to -1.0 are UM 
420, UM 448 and Mkn 1089. 

In Table 121 we present parameters of the H ii region on 
the emission extracted with the aperture l'.'5x3'.'0 on the 
MMT spectrum, and with the aperture 2"x2'.'8 - on the 
BTA spectra. They include both the observed line inten- 
sities and intensities corrected for interstellar extinction, 
relative to that of H/3 along with the equivalent width 
EW(H/3), observed flux of the H/3 emission line, and ex- 
tinction coefficient C(H/3). To correct for extinction we 
used the Whitford fl9585 reddening law for our Galaxy. 

The derived values of Tc(Oiii), Te(Oii), rc(Siii) and 
Nc(Sii) for this supergiant Hii region are shown in 
Table 121 The ionic and total abundances in the Table are 
derived from the two MMT spectra according to the meth- 
ods described in ITL94 and in Izotov et al. H1997I here- 
after ITL97). Since they showed a mutual consistency, we 
averaged their derived parameters. 

In the derivation of the O, N and Fe abundances, 
the broad low intensity underlying components of the 
strongest emission lines Ha, H/? and [O ill] AA4959, 5007, 
and WR broad lines were taken into account (see Sect. 
I3.2|l . For the abundance determination only the fluxes 
of the narrow (instrumental width) components of these 
strong lines were used. In the derivation of the CI abun- 
dance, we used the respective formula from AUer (1984) 
with the ICF adopted from Stasiiiska (,1990,1 . The abun- 
dances found here are well consistent with the previous re- 
sults from BTA those in Kniazev et al. H2000a|l and those 
based on the SDSS spectra (Kniazev et al. I2003|l . 

As evident from Table O the object belongs to the 
group of the most metal-poor H ii-galaxies, which in- 
cludes only 1-2 per cent of the known BCGs. Notice the 
very high abundance of nitrogen. The N/0 ratio is ^^6 
times higher than the typical value for all known XMD 
BCGs. As shown in IT99, BCGs with 12+log(O/H)<7.60 



have log(N/0) «-1.60 with a very small scatter (Fig. 
01). The relative abundances of the other elements in HS 
0837+4717 are close to the mean values for the XMD 
BCGs (IT99). 

3.2. WR features and broad components of strong 
emission lines 

The Wolf-Rayet features of the blue bump (near -4650 A) 
and broad low-contrast components of the strong emission 
lines H/5, Ha, [O ill] AA4959, 5007 are evident in both the 
MMT and BTA spectra of HS 0837+4717. This allows us 
to check the accuracy of the derived parameters. 

3.2.1. WR blue bump features 

In Figure 121 we show an expanded view of the BTA spec- 
trum with the Gaussian decomposition of the narrow lines 
and the broad component of He ii A4686 and the WR lines 
Niii A4640, attributed to late WN stars, and Civ A4658, 
attributed to WC stars (e.g., SV98). The latter also could 
appear as the red bump centered at A5808. This bump, 
however, is barely detected in our spectra. Its flux as 
estimated (S/N ratio ^^2) from the MMT spectrum, is 
2.4x10""'^^ erg cm~^ s~^. The parameters of fitted lines in 
different spectra are consistent within their rather large 
uncertainties. The parameters of narrow and broad lines 
are presented in Tables (21 and 0] respectively. 

3.2.2. The broad components of the strong emission 
lines 

The broad components of the strong emission lines of Ha, 
H/3 and [Oiii] AA4959,5007 originate in the ionized gas 
moving with velocities of ^^1-2 thousand km s~^. Such 
high velocities are related either to supernovae shells or 
to stellar winds of massive stars. Since the widths of 
the Balmer and the [O ill] lines are similar, it is reason- 
able to suggest that both originate in the same regions. 
While broad features centered at H/3 (the He ii A4861 +H/3 
blend) and Ha (the Hen A6562 + Ha blend), charac- 
teristic of WR stars (e.g., SV98), should contribute to 
the broad components of the Balmer lines, their rela- 
tive strengths calculated from models are expected to be 
tens of times lower than ones observed in HS 0837+4717. 
Therefore their origin in this galaxy, as well as in other 
similar cases, can probably be attributed to rather young 
remnants of multiple SN explosions. Having the line inten- 
sity ratios for the broad components, we can address the 
physical conditions in the regions of their formation. The 
parameters of broad components of these four lines are 
presented in Table 01 It is evident that these components 
have systematically higher line fluxes in the BTA spectra. 
This can be attributed to the wider slit (2" at BTA vs 1'.'5 
at MMT) and to the different PA. 

The Balmer decrement of the broad line components, 
/(Ha)//(H/3) appears unusually large. Accounting for the 
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significant uncertainties, we accept the average of the 
MMT and BTA ratios for the fluxes of the broad and nar- 
row components of H/3, equal to (2.7±0.5)%. The same 
average ratio for Ha is (19±1.5)%. Adopting a flux ratio 
of the narrow Ha to H/3 lines, corrected for extinction, of 
2.76, a Balmer ratio for the broad components is 19.4±4. 

Such a large Balmer decrement is rather unusual. To 
ensure its reality we checked two possible factors, which 
could in principle enhance the Balmer ratio. The first one 
is underlying stellar absorption in the Balmer lines, which 
has an EW comparable to the EW(broad H/3), and hence 
could affect it. However, since the FWHM of stellar ab- 
sorption lines (^200 km s~^) is typically ^7 times smaller 
than that of broad emission component, the effect of stel- 
lar H/3 absorption will decrease the flux of broad emission 
H/3 at most by 15%. The second factor could be purely 
instrumental: light scattering in the spectrograph of the 
strong lines could produce some artifacts with arbitrary 
ratio of intensities. This certainly can be excluded, since 
we observed another object during the same night at the 
MMT and found even stronger emission lines, but with no 
hints of the broad components. 

Thus, this very large Balmer decrement is either re- 
lated to substantial internal extinction of the regions 
where the broad Balmer lines are formed, or to the ex- 
citation of hydrogen levels by electron collisions, or both. 
For the latter, Ha/H/3 ratio can reach the value of 5.8 
for Te^lOOOO K (Chamberlain HHSSIl- Thus, even for the 
case of electron excitation, some significant extinction (ac- 
counting for the additional factor of at least ~3.3±0.7 
in the Balmer line ratio) should be invoked for the re- 
gions of the broad Balmer lines formation. Using the 
Whitford extinction law, the latter possibility corresponds 
to E{B — V)—1.16zt0. 2 and an additional extinction in the 
B-band of Ab~5±1 magnitudes! The latter indicates that 
the broad line emission forms in, or propagates through, 
a very dusty medium, contrary to what is observed for 
the narrow line Hii region type emission. If no alterna- 
tive explanation of the large Balmer decrement is found, 
this implies that the fluxes of broad features, and of the 
unknown optical continuum related to them, should be 
corrected by a factor of 25-100, depending on the wave- 
length. Such correction will raise the total luminosity of 
the BCG by a factor of several. The highly obscured re- 
gions are quite typical of starbursts observed in mergers. 
Moreover, in the XMD BCG SBS 0335-052 E mid-IR ob- 
servations reveal a highly obscured starburst (Hunt et al. 

mm . 

We also point out the very large difference between 
the ratio of /(A4959)//(H/3) for the narrow and the broad 
components: for the average of the MMT and BTA data 
they are 1.88 and 0.62, respectively. Is this due to the effect 
of significant coUisional de-excitation of 0++ ions, or to a 
significantly lower O /H ratio in the regions of broad line 
formation, or their lower Tc? It is possible that all these 
factors combine to produce the observed effect. To answer 
these questions, higher S/N ratio data will be helpful. 



3.3. Spatial structure of the emission regions 

Below we exploit the advantage of long-slit spectra to ex- 
amine various morphological features of the object. The 
seeing during the MMT observations {9=1'.'6) and BTA 
observations in the red {9=1'! 3) was significantly less than 
the full extent of both line and continuum emission (8"- 
10"), and allows the performance of a coarse analysis of 
the BCG substructure. Note that what we refer to as 
seeing is the width of an unresolved object as measured 
perpendicular to the spectrograph's dispersion and so in- 
cludes contributions from the atmosphere as well as the 
telescope and spectrograph optics. 

In Figure[Slwe show the intensity profiles of the contin- 
uum along the slit on the MMT spectrum averaged over 
the line-free spectral ranges (in the rest frame system) 
of 4110-4310 A and 5900-6400 A. For comparison, the 
intensity distribution in a narrow band, centered on the 
Ha+[Nii] lines is also presented. 

The intensity distributions of the continuum in Fig. El 
and of the strong emission lines are different. The latter 
have weU fitted PSF profiles, with the FWHM close to 
that of seeing. Some very faint wings/pedestal are present, 
which have just the appearance of the PSF at a few per 
cent level of the peak value. 

The spatial distribution of continuum in the BCG 
shows some additional features. The main peak, related 
to the current starburst, fully coincides in position (along 
the slit) with the peak of the line emission. This starburst 
continuum component dominates the other components. 
However, an additional component (knot) on the MMT 
spectrum is seen at '^2'.'2 (1.8 kpc) to the SSW, with the 
relative fiux of ~25% in the blue and of ~28% in the red. 
This imphes that the {B - R) of the faint 'knot' is -0™13 
redder than that of the current starburst, which domi- 
nates the BCG light and determines its integrated colours. 
Based on the relation above we can approximately esti- 
mate the {B — R) colour of the faint knot. From the data 
in Table [SJ the integrated (_B — i?)o is 0.15. This colour is 
affected by the strong emission lines included in both fil- 
ters. Accounting for the EWs of these lines and for the ef- 
fective widths of both filters, we estimate that the (B — R) 
of the continuum component is redder due to the emission 
lines by ~0.12. This implies that the net continuum of HS 
0837+4717 has a {B-R)o '-O.OS. The latter is very similar 
to the colour {B — R)=-0.05 of the mixture of a starburst 
with an age of 3.7 Myr and the continuum produced by gas 
at To=20000 K. The former, according to the PEGASE.2 
model (with a Salpeter IMF with limits of 0.1 and 120 
Mq and Z=Zq/20), is (B - i?)=-0.24, while the latter, 
from the data in AUer (19M), has {B - R)=+0.53. They 
are mixed in the proportion that results in the observed 
narrow-line EW(Ha)~1400 A. We accept the bright com- 
ponent continuum colour of (B — i?)=-0.05, and from the 
latter we obtain for the faint knot (B - i?)[f™* ^0.08. For 
an instantaneous star formation (SF) with the same IMF 
and Z=Z0/2O, this (B — R) corresponds to an age of ~25 
Myr. Apart from this fainter knot of continuum emission, 
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there is also some diffuse component, stretching further to 
the SW and NE, as seen at the F-band image in Fig. ^ 

We estimate the intrinsic size of the supergiant H ii re- 
gion from a comparison of the FWHM of its strong emis- 
sion lines with that of PSF. Unfortunately, for the sharpest 
images {d—l'!3, BTA observations in the red), the seeing 
fluctuated, therefore it is difficult to accept the value of 
the FWHM for the PSF. For the MMT observations the 
situation is better. A nearby star appears on the slit and 
can be used to estimate the PSF at various wavelengths. 
The weighted mean of the FWHM(PSF) in red and green 
is 1'.'60. The FWHM (intrinsic), derived as the weighted 
mean of the values in red and green, is ^^O'.'TitO'.'O?. 
FWHM(intrinsic) were calculated by quadratic subtrac- 
tion of the FWHM(PSF) from the measured FWHM of 
strong lines. The size of the faint knot is difficult to mea- 
sure directly. However, from the decomposition of the light 
distribution along the slit, shown in Figure the FWHM 
of this knot is ~2'.'2. Adopting the same FWHM of the 
PSF, we obtain a coarse estimate for the intrinsic size 
of the faint knot as FWHM (knot )~1'.' 5. These sizes corre- 
spond to a FWHM linear diameters of the two components 
of ~0.58 kpc and ^^1.25 kpc, respectively. 

Summarizing this analysis, we find that along the ma- 
jor axis at its brightest part, this BCG has a double- 
nucleus structure with characteristic sizes of the compo- 
nents of 0.6 and 1.25 kpc, respectively, and a projected 
distance of 2'.'2 (1.8 kpc) between the components. The to- 
tal extent of the underlying low-surface-brightness (LSB) 
component is >10" (8 kpc). The nature of the faint knot 
is not clear. But since it is also quite blue, one can expect 
that this is also related to the recent SF event. 



3.4. Ionized gas kinematics 

We performed an analysis of the gas kinematics follow- 
ing the method described by Zasov et al. I|2()[)()|l . The 
Ha line position-velocity (P-V) diagrams along the slit 
at PA=44° and 130° are shown in Fig. |H1 The first one is 
close to the major axis of the outermost isophotes, while 
the second is nearly perpendicular. The full extent of Ha 
emission as seen from these spectra is the same in both 
directions: 9"-10" (or ~7-8 kpc). The appearance of P-V 
diagrams, however, is very different. For PA=130° the ra- 
dial velocity is practically constant for a range from -3" 
to ~3'.'5, or about 65% of the entire extent of the observed 
Ha-emission. This is consistent with a negligible expected 
projection of the rotation velocity for this PA. The veloc- 
ity curve at the galaxy edges looks, however, somewhat 
disturbed. This could be a result of a recent interaction or 
some indication for a gas blowout, confined in some direc- 
tions. Better S/N ratio data for these regions are necessary 
to confirm the significance of the velocity disturbances. 

The slit position at PA=44° is expected to be the opti- 
mal one to see the overall rotation of this galaxy. However, 
the respective P-V diagram does not look like a typical 



Table 3. Abundances in the supergiant H ii region 



Value 



MMT 



BTA 



re(oiii)(K) 

To (Oil) (K) 

re(siii)(K) 

N,{SIl){cm-^) 

0+ /U+ (xlQ-^) 

0++/H+(xlO-5) 

O+++/H+(xl0-5) 

0/H(xlO-^) 

12+log(0/H) 

N+/H+(xlO~'') 

ICF(N) 

log(N/0) 

Ne++/H+(xlO-5) 

ICF(Ne) 

log(Ne/0) 

S+/H+(xlO-^) 
S++/H+(xlO-^) 
ICF(S) 
log(S/0) 

Cl++/H+(xlO-^) 

ICF(Cl) 

log(Cl/0) 

Ar++/H+(xlO-^) 
Ar+++/H+(xlO~^) 
ICF(Ar) 
log(Ar/0) 

Fe++/H+(xlO-^) 
ICF(Fe) 
log(Fe/0) 
[O/Fe] 



18,500±700 
15,200±500 
17,000±500 
300±300 

0.345±0.035 
3.941±0.355 
0.082±0.014 
4.369±0.357 
7.64±0.04 

5.157±0.701 

12.65 
-0.83±0.07 

0.610±0.058 

1.11 
-0.81±0.05 

0.695±0.061 
3.437±1.294 

2.91 
-1.56±0.14 

0.121±0.004 

2.19 
-3.22±0.14 

1.003±0.101 
1.402±0.193 

1.01 
-2.26±0.05 

0.663±0.321 

15.79 
-1.62±0.21 
0.20±0.21 



18,300±800 
15,100±600 
16,900±600 
400±300 

0.357±0.042 
3.969±0.419 
0.095±0.022 
4.421±0.422 
7.65±0.04 

5.265±0.660 

12.40 
-0.83±0.07 

0.606±0.068 

1.11 
-0.82±0.06 

0.767±0.074 
5.196±1.282 

2.87 
-1.41±0.10 



1.395±0.317 

1.71 
-2.27±0.11 

0.622±0.212 

15.45 
-1.66±0.15 
0.24±0.15 



rotation. We discuss below two possible explanations for 
the unusual distribution of the radial velocity. 

The first one treats the P-V diagram as the conse- 
quence of counter-rotation. At positive coordinates along 
the slit we see the receding branch of the rotation curve 
with a velocity amplitude of 50-70 km s~^ (depending on 
how the real curve is drawn between the two last points) 
relative to the brightness peak velocity. However, the sym- 
metric (approaching) continuation of this branch at nega- 
tive coordinates is not seen. Instead, the receding branch 
is present, with an amplitude of 20-30 km s~^. This could 
probably indicate 2 dynamically decoupled gas systems. A 
difference of ~30 km s^"'^ between the radial velocities near 
the peak of the Ha emission for PA=130° and PA=44° can 
be attributed to small difference in the zero-points of the 
dispersion curves for the two independent spectra. 

Another explanation treats the part of P-V diagram 
at negative coordinates as a significant velocity distur- 
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Fig. 5. Left panel: Intensity distribution along the slit of the continuum between Hi5 and H7. The dashed line shows the 
PSF (point-spread function) fitted to the central peak (related to the starburst). The dotted line shows the difference 
between the observed full light and that of the central bright compact region. The additional component, centered at 
Ri2'.'2 to the SSW is ^^0.25 of the bright peak amplitude. Middle panel: Same as in the left panel, but for continuum 
in the red. The relative flux of the fainter knot is 0.28, that is a factor of 1.13 larger than in the blue. This implies 
that this knot in {B — i?) is ~ 0™13 redder than the central starburst, which mainly defines the integrated colours of 
the BCG. Right panel: Same as the left panel, but for intensity centered on the Ha line. No significant signal is seen 
at the periphery besides the scattered light from the central bright region. 



bance relative to the extrapolated approaching branch. 
This could be related to a strong outflow from the recent 
SF episode. If one is to draw the approaching branch of 
the rotation curve as a symmetric extrapolation of the re- 
ceding one, the resulting difference between the observed 
data and drawn approaching branch shows us the form 
of the disturbance. This resembles the common wave-like 
form of supershells, as they are usually seen on P-V dia- 
grams of starbursting galaxies. The velocity of such super- 
shell relative to the background rotation velocity is ~70- 
100 km s~^. Its visible extent is ^2". The results could 
indicate that only a half of the shell is seen, which is closer 
to the center. Therefore, its full extent could be ~4", or 
'^B.B kpc. Such features in P~V diagrams are detected in 
many BCGs, and the most plausible explanation is related 
to the expansion of supershells formed in the ionized gas 
due to the recent star-forming activity in these objects. 
Supershells are produced by hot bubbles, originating from 
the energy injection of numerous massive star winds and 
supernovae (SNe) explosions into the galactic ISM (e.g., 
Tenorio-Tagle & Bodenheimer 1198^ . 

In particular, a similar, but scaled down according 
to the parameters of the host galaxy, supershell in HS 
0822+3542 was discussed by Pustilnik et al. pK)3afl . The 
energy of that supershell with a diameter of 0.48 kpc and 
a velocity amplitude of 30 km s~^ was provided by 13 
SNe. Assuming the same density of ambient gas of 0.1 
cm~^, with the supershell parameters as observed in HS 
0837+4717 (adopting a minimal T4hoU=70 km s~^), we de- 
rive a total energy larger by a factor of ^^1770. This cor- 
responds to ~23000 SNe exploded during a period of the 
order of one dynamical age of the supershell. The latter 
is idyn = 0.6 i?shcii/V^shcii=13.7 Myr (Weaver et al. [W77|l . 
This would imply the existence of a star cluster with an 
age of ~17 Myr, provided the mean rate ~1700 SNe per 



Myr during last ~14 Myr (recall, that SNe begin to ex- 
plode in the cluster after ~ 3 Myr). From this SN rate 
and a PEGASE.2 model with the same IMF as above, we 
derive the mass of such a cluster of ^7.6x10*^ Mq. 

For a galaxy of this high luminosity and an appar- 
ent diameter of ^10 kpc, the estimated rotation velocity 
is rather small. The inclination correction, if applicable, 
is modest. From the apparent axial ratio of this BCG 
{p — h/a = 0.67), it follows, that in the case of a stable 
rotating disk, i is quite large. Indeed, from the well known 
formula: cos^(i) = (p^ — q^)/ (1 — q^), and the range of in- 
trinsic ratio for disk-like galaxies ofq = b/a = 0.2—0.3, the 
i would fall in the range of 49°-52°. The respective correc- 
tion of l/sin(i)=1.27-1.32. In this case the corrected I^ot 
can reach the value of 65 to 90 km s~^, within the range 
characteristic of low-mass galaxies. 

If this object is a merger, gas motions due to dissipative 
collisions could be significantly decoupled from the motion 
of stars. Until the system relaxes, the gas motions are bad 
tracers of mass distribution (e.g., Amram & Ostlin l2()()l|l . 
Therefore the amplitude of the rotation velocity observed 
in Ha could be unrelated to the overall gravitational po- 
tential. 

3.5. Imaging and photometry 

In the light of the discussion on the possible merger nature 
of HS 0837+4717, it is important to inspect the 2D mor- 
phology of the galaxy. Its grey-scale T/-band image with 
isophotes superimposed is shown in Figure Q 

The overall structure is asymmetric and somewhat dis- 
turbed, both near the center, and in the peripheral regions. 
In particular, the fainter component at ~2" to the South 
(which is seen in the continuum on the long-slit spec- 
trum) is apparent on this image, however smeared due to 
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Fig. 6. Left panel: Intensity distribution of Ha along the slit at PA=44° and the respective P-V diagram (binned with 3 
pixels, or 1'.'2, corresponding to the seeing). The radial velocity increases monotonously in the SW direction, resembling 
the receding branch of rotation curve with the amplitude of 50-70 km s^^. However, the expected approaching branch 
for negative coordinates in the direction to NE is absent (see text for discussion). Right panel: Same as in the left 
panel but at PA=130°. The velocity variations in the central part of this cross-section are significantly smaller than 
those at PA=44°. This would be expected due to the small projection of the rotation velocity onto the minor axis. 
There seems to be a large velocity disturbance at the NW edge, where radial velocity steeply raises to 60-70 (±30) 
km s~^ in the two last points. The small velocity variations are also observed on the opposite side of the galaxy. 



the seeing effect. The outermost isophotes indicate several 
appendages to the main body. Disregarding the faintest 
isophote as being subject to noise, we notice two small 
features, at NNW and SSE of the center, the latter being 
much more prominent. The BCG images available from 
the Data Release 1 (DRl) database of the Sloan Digital 
Sky S urvey (SDSS) (York et al. 2000, Abazajian et al. 
I2003|l corroborate this disturbed morphology. 

As described above, small transparency variations took 
place during imaging with the Wise 1-m telescope. In or- 
der to minimize their effect on the object's integrated mag- 
nitudes, we selected from each of the nights the bright- 
est value in each of the bands, and compared them with 
the magnitudes, derived through the convolution of the 
MMT spectrum. For [/-band we could not use an MMT 
measurement due to the limited wavelength coverage. 
Therefore, even though the [/-band magnitudes for both 
nights are consistent within the observational uncertain- 
ties, they should be treated with caution. The brightest 
values for the Wise measurements in B and R are consis- 
tent with the MMT data within the uncertainties, there- 
fore we accept them with an r.m.s. error of 0™10. In the 
y-band, however, the brightest magnitude on Wise data 
is still 0V^27 fainter than the respective magnitude derived 
from the MMT spectrum. Therefore, for V we accept the 
MMT derived value. The integrated photometry data are 
summarized in Table |S1 

The very red {B~V) (-0.9) and blue (V - R) ( — 0.7) 
colours of this object are related to the very large EWs of 
the [O III] A4959 and A5007 A emission lines, redshifted to 



the center of l/-band. This can be checked by the convo- 
lution of the deredshifted BCG spectrum (see Table 0. 

3.6. Hi results 

Effectively, we observed HS 0837-1-4717 in "on-l-ofF' posi- 
tions for 320 minutes and, after smoothing, obtained an 
r.m.s. noise of 1.3 mK. This translates to an r.m.s. in flux 
density of 1.6 mJy. Since no signal was detected, we put 
a 2(7 upper limit of 3.2 mJy on the peak flux density from 
the galaxy. The full amplitude of ionized gas velocity vari- 
ations, ~100-140 km s~^ is used to estimate the upper 
limit of the H i flux. Adopting this value as an estimate of 
14^50 for the Hi profile, the upper limit of the integrated 
Hi flux is of (0.32-0.45) Jy km s^^. The respective upper 
limit on theHimassisM(Hi,BCG) < (2.1-3.0) x 10^ Mq. 

4. Discussion and conclusions 

4.1. General comments 

The integrated parameters of HS 0837-1-4717 are summa- 
rized in Table |S1 There are two unusual properties, which 
are probably related. The flrst one is the high luminosity 
of this BCG for its very low metallicity (see, e.g., plots of 
0/H v ersus Mb in Kimth & Ostlin 121001 Pustilnik et al. 
I2003b|l . The second one is the very large nitrogen abun- 
dance excess: N/0 is a factor of '--^6 larger than for other 
BCGs with similar 0/H. The flrst property can be at- 
tributed to an exceptionally strong starburst, which causes 
an additional brightening by several magnitudes. Pustilnik 
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Table 4. Parameters of the WR lines and of the broad 
components 



Table 5. Main parameters of HS 0837+4717 





MMT 


BTA 


Flux (10""" ergs s^'cm"^)" 




4640 N III (WR) 


2.5 


3.8 


4658 C IV (WR) 


2.6 


5.0 


4686 He ii (WR) 


2.4 


3.7 


4861 H/3 


6.3 


15.0 


4959 [0 III] 


11.6 


18.3 


5007 [0 III] 


34.8 


55.0 


6563 Hq 


159.4 


261.3 




Equivalent width [A] 




4640 N III (WR) 


1.4 


2.1 


4658 C IV (WR) 


1.5 


2.8 


4686 He ii (WR) 


1.4 


2.1 


4861 H/3 


4.0 


7.8 


4959 [O III] 


7.9 


9.8 


5007 [0 III] 


24.3 


30.0 


6563 Hq 


238.4 


302.2 




FWHM [km s~"] 




4640 N III (WR) 


1300 


1300 


4658 C IV (WR) 


1300 


1300 


4686 He ii (WR) 


1300 


1300 


4861 H/3 


1500 


1500 


4959 [O III] 


1500 


1500 


5007 [0 III] 


1500 


1500 


6563 Hq 


1300 


900 


Broad/Nebular flux ratio [%] 




4686 He ii (WR) 


36 


46 


4861 H/3 


1.9 


3.5 


4959 [0 III] 


1.7 


2.2 


5007 [0 III] 


1.7 


2.2 


6563 Hq 


16 


22 




WR/O stars ratio 




Observed 


0.02 




Model*" 


0.02 





" Flux errors are ~50% for the broad He ll A4686 line and 
~(15-25)% for the broad components of the strong lines. 
'' data from Schaerer & Vacca I1998II . 



et al. I|2001b|l argued that the majority of starbursts in 
BCGs are triggered by interactions of various strengths 
with other galaxies. The range spreads from weak tidal 
interactions through strong tidals and the sinking of satel- 
lites to merging of low-mass galaxies. The latter case can 
be considered as a down-scaled analog of ultraluminous 
IR galaxies, mergers of massive gas- rich galaxies (e.g., 
Genzel et al. I2()()l|l . Below we discuss the possibility that 
HS 0837+4717 as well could be an advanced merger. 

4.2. Possible evidences of merger 

Tellcs and Terlevich l)1997|l found that the most lumi- 
nous BCGs show disturbed external morphologies, while 
low luminosity BCGs show regular morphology. In their 
detailed studies of several luminous BCGs, Ostlin et al. 
H1999ll^ni7T|l found morphological and kinematic evidences 



Parameter 



Value 



R.A.(J2000.0) 
DEC.(J2000.0) 

^tot 

(u - i3)r 

{B - V)r 

(v - i?)r 

VhoI (km s"^) 
Distance (Mpc) 

Opt. size (")* 
Opt. size (kpc) 
12+log(0/H) 
re(Oiii) (K) 
iVe (cm-3) 
log(N/0) 
Hi int.flux(5) 
M(Hl) (10^ Mq) 
M(Hi)/Lb'^' 



08 40 29.84 
+47 01 09.1 

0.11 
18.18+0. 10'^> 
-0.41+0.14'^' 

0.86+0. 14'^'^' 
-0.71+0.14^"'^^ 
12630<^> 
170.5(2) 
-18.1(2) 
~8x7(2) 

6.5x6(2) 
7.64+0.03(2) 

18000(2) 

<400(2) 
-0.83+0.05(2) 
<0. 31-0.43(2) 

<2.5(2) 

<0.9(2) 



tot corrected for redshift is 0^21. 
tot corrected for redshift is -0™02. 



(1) - (B - l/)i'< 
{l)-iV-R)l 

(2) - this paper. 

(3) - corrected only for the Galactic extinction, 

(4) - a X fe at fiB =25 mag/n". 

(5) - for Wso <100-140 km s-\ in Jykm s"^ 

(6) - in units (M/Lb)© 
(N) - data from the NED. 



for the merger nature of these galaxies. HS 0837+4717 is 
similar to their BCGs in several aspects. 

As shown above, the overall disturbed morphology of 
this BCG, its " double- nucleus" appearance, and the pres- 
ence of appendages on the periphery, presumably imply 
some recent disturbance. However, no disturbing galaxies 
are visible in the environment of this BCG (see Sect. l4.5|l . 

In addition, the distribution of the ionized gas ve- 
locity also appears disturbed, mainly near the edges of 
the galaxy. The gas motions along the apparent major 
axis show large deviations from the expected regular ro- 
tation. These could be related either to a counter-rotating 
gaseous component, or to a very energetic supershell. The 
two knots, related to the current/recent starbursts are 
very massive. Indeed, from the photometry in table |S1 and 
numbers derived in Sect. 13.31 from a PEGASE.2 model 
with Z=Zq/20, and a Salpeter IMF as above, we find 
the following. For the brighter component with Mb=- 
17.9 and an age of 3.7 Myr the mass of starburst is 
0.8x10^ Mq. As shown above, a significant part of this 
starburst could be highly obscured by dust. Therefore, its 
unobscured blue luminosity could be several times higher 
than that listed here. This would translate to a total mass 
for this starburst region of >3xl0^ Mq. For the faint 
knot A'/B,faint=-16.4. Its colour {B — R) corresponds to 
an instantaneous SF burst with an age of 25 Myr. With 
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the same IMF as above, PEGASE.2 predicts a star clus- 
ter mass of 2.7x10^ -^^o- Also, from the parameters of 
the suggested supershell at the NE edge of the galaxy, we 
estimate the mass of the embedded star cluster with an 
age of '^IT Myr as '--^0.76x10'' Mq. Such massive star- 
bursts should involve molecular clouds with the masses of 
10^ Mq or more. We will return to this point in Sect. 14.41 
Finally, as shown, e.g., by Mihos & Hernquist (|296), very 
energetic starbursts are characteristic of merging galaxies. 

An alternative option for such a case is a strong in- 
teraction with an undetected massive object, such as very 
LSB galaxy or a hydrogen cloud devoid of stars, similar 
to the case of Dw 12254-0152 (Salzer et al. [T!M|l . 

4.3. WR star population 

The WR stellar population is observed only in four of 
about 30 XMD BCGs known. They include I Zw 18 
(Izotov et al.^^ Legrand et a,]. [T#7|l . SBS 0335-052E 
(Izotov et al. lMTT)! . Tol 1214-277 (Fricke et al. 2001) and 
UM 133 (Kniazev et al.EOOl, 2004, in prep.). Thus, HS 
0837-^4717 is the fifth XMD BCG with a detected WR 
population. While the accuracy of the flux measurements 
of the broad WR features on the available spectra is far 
from sufficient for accurate prediction, it is useful to ob- 
tain an estimate of the WR-to-0 star relative numbers. 

To derive iV(WR)/iV(0-hWR), we foll ow the approach 
described in detail by Guseva et al. 12000(1 . The number of 
WR stars was derived from the luminosity of the WR lines. 
Their fluxes were derived directly from the 1-D spectra, 
as discussed in Sect. 13.2.11 No aperture correction was 
applied to the WR-line fluxes, assuming that this emission 
is produced in a compact region with an angular size less 
than the width of the slit. The number of O stars was 
derived from the total H/3 luminosity for the brightest SF 
region (icor), corrected for the contribution of WR stars 
and for the aperture effect. For the latter we adopt the 
circular symmetry of the bright SF region as seen in the 
direct image. To estimate the fraction of the undetected 
H/3 flux, the radial distribution of the H/3 line intensity, 
shown in Figure El was used. The respective correction 
factor for the total H/3 flux is 1.66. 

The number of O stars is derived from the number of 
ionizing photons Q™"^, which is related to the total lumi- 
nosity of the H/3 emission line icor(H/3) by 



assumed that the average Lyman continuum flux per one 
WR star is comparable to Qq'''^ and is equal to 1.0 

X 10^^ s-i (Schaerer et al. fTM?^ . Finally 



,(H/3) =4.76 X 1{)-^^QI° 



(3) 



For a representative 07V star we adopt the number of 
Lyman continuum photons emitted to be Qq''^ = 1 x 10"*^ 
s~^ fLeitherer ll990|l . The total number of O stars is then 
derived from the number of 07V stars by correcting for 
other O star subtypes, using the parameter 770 introduced 
by Vacca & Conti H1992|l and Vacca fl994|. 770 was calcu- 
lated in SV98 as a function of the time elapsed since the 
onset of an instantaneous star formation burst, as inferred 
from the EW{B.(3). The Salpeter IMF with the Afiow and 
Mup limits of 0.8 Mq and 120 Mq was adopted here. We 



N{0) 



(4) 



In Sect. 13.2.11 we performed the estimates of the ob- 
served fluxes of all broad WR lines in the blue and red 
bumps (see Table 0J). The latter are used to get the lu- 
minosities in the respective lines. Then, the luminosity of 
red bump is used to estimate the number of representative 
WC4 stars, having ^5808=3. Ox lO^^ erg s^^: AVc=280. 
The number of representative WN7 stars can be estimated 
either from its adopted luminosity in broad He II line, or 
from the sum luminosity of the blue bump, subtracting the 
contribution of CivA4658, caused by WC stars. The un- 
certainties of all WR broad features are large, about 50%. 
Therefore the derived numbers of WR stars are only in- 
dicative. Better signal-to-noise spectra are needed to study 
the WR population in this BCG at the quantitative level. 

Assuming the luminosity of a single WN7 star to 
be 1.6 xlO'^^ ergs s~^ in the broad He 11 A4686 line or 
1.99x10'^^ ergs s~^ in the blue bump (SV98), we estimate 
from the WR line fluxes derived on the MMT spectrum, 
iVwN to be either 500, or 860. Adopting the average of the 
two estimates for WN stars, we get the sum number of WR 
stars (WN-I-WC) in the young starburst to be '-^1000. 

Following to SV98 and using the observed value of 
EW{Il(3) = 230 A, we estimate the burst age and the pa- 
rameter rjo (updated model version of 1999, as presented 
in Starburst99) in HS 0837+4717 as 3.7 Myr, and -0.65, 
respectively. The total number of O stars in HS 0837+4717 
is found to be N(0) w 47000. This yields a relative num- 
ber of WR stars to O stars iV(WR)/(iV(0 + WR) of 0.02, 
coinciding with the model prediction in SV98. 

4.4. A large excess of the nitrogen abundance 

The possibility of enhanced N abundance in galaxies with 
the spectral signatures of WR stars, in comparison to 
other galaxies of similar metallicity was first noticed by 
Pagel et al. H1986(l where no details on the strength of this 
effect were given. The occurrences of a significant nitro- 
gen excess in galaxies are very rare. For example, among 
more than one hundred BCGs with well measured heavy 
element abundances (e.g., IT99; Ugryumov et al. .200311 
only four galaxies (Mkn 996, UM 420, Mkn 1089 and UM 
448) exhibit a nitrogen excess of a factor of 3 and larger, 
relative to the value appropriate for their 0/H (Thuan et 
al. nMHI IT99). All four have 12+log(0/H)>7.9 and show 
broad WR lines (Guseva et al. 2000 Thuan et al. [M?H| . 
One more similar galaxy, Haro 11 is found among the lu- 
minous BCGs studied by Bergvall & Osthn With 
12+log(0/H)=7.9 it has a nitrogen abundance excess of 
a factor of —6. Strong signatures of WR population are 
detected in this BCG as well. 

An excess of the nitrogen abundance by a factor of 
3 is also detected in two central Hii regions (both with 
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Table 6. Star burst galaxies with significant nitrogen excess 



lAU Name 


Synonyms 




Dtot 


Mb 


12+log(0/H) 


Alog(N/0) 


Notes 


0034-3349 


Haroll 


6175 


14.6 


-20.0 


7.90 


0.8 


merger 


0125-061 


MKN996 


1622 


15.1 


-17.1 


8.00 


0.6-1.4 


probable merger 


0218+003 


UM420 


17514 


16.5 


-20.3 


7.89 


0.5 


sinking merger? 


0459-043 


MKN1089 


4107 


15.0(13.3) 


-19.1(-20.8) 


8.07 


0.5 


bright knot in interacting gal. 


0837+4717 


HS,PC 


12630 


18.2 


-18.1 


7.64 


0.77 


merger? 


1139+006 


UM448 


5560 


14.7 


-19.9 


7.98 


0.5 


2 knots in center - merger? 


1337-313 


NGC5253 


440 


10.9 


-17.4 


8.16 


0.5 


Im pec, starburst 



° data from the NED. 



detectable WR population) of the nearby starbursting 
dwarf irregular galaxy NGC 5253 with 12+log(0/H)=8.16 
(Kobulnicky et al. 119971 and references therein) . In Table 
El we summarize the parameters of all known such galax- 
ies, and below draw some tentative empirical correlations. 
First, most of these galaxies could be identified as various 
stages of mergers. Second, as a subgroup of starbursting 
galaxies, 70% of them are luminous (Mg <-18™0), thus 
are not classical dwarf galaxies. Third, all these objects 
are WR galaxies. Fourth, all these BCGs but one have 
12+log(0/H) in the range of -7.9 to 8.16. HS 0837+4717 
is the only such galaxy belonging to the XMD group. 

The scenarios of localized nitrogen enrichment, as well 
as the constraints imposed by the known mechanisms of 
fresh elements mixing and dispersal were discussed in de- 
tail by Kobulnicky et al. ( ISSTJ (see also Roy & Kunth 
I1995|l . They concluded that, in the light of the expected 
short timescales for dispersal in the interstellar medium 
(ISM) of a few Myr, the localized nature of the N en- 
richment suggests a very recent pollution event, proba- 
bly connected with the onset of massive star winds. The 
alternative analysis of Tenorio-Taglc (1996) suggests the 
timescale of —100 Myr to mix fresh metals from a cluster 
of massive stars with the gas in a disk galaxy. However, 
this can be barely applicable for very complex flow pat- 
terns expected in the vicinity of merger starbursts. In par- 
ticular, the very rare occurrence of objects with a large N 
overabundance, similar to that observed in Mkn 996 and 
HS 0837+4717, favors the general idea of the both short- 
time scales: for localized N pollution and its fast dispersal. 

It is worth noting that a selection effect, related to 
the starburst strength and its trigger mechanism, will af- 
fect the observed frequency of this phenomenon, probably 
determining the overall scale of the polluted region. In 
particular, the large N excess in the central starbursts of 
NGC 5253 would be barely detected with ground-based 
spectroscopy, if this galaxy was a few times more distant. 

4.5. The global environment 

On the large scales, the volume around HS 0837+4717 
is well sampled only by the objects from the Updated 
Zwicky Catalog (Falco et al. I1999|l , which is complete to 
B =15™5. At the BCG distance they all are luminous. 



with Mb <-20™6, and are located outside the sphere with 
R = 8.4 Mpc, centered at HS 0837+4717. This implies that 
HS 0837+4717 is situated in a void (as several other XMD 
BCGs are, Pustilnik et al. I2003ajl . New data on the dis- 
tribution of fainter galaxies with measured rcdshifts from 
the SDSS DRl database (Abazajian et al. .200311 around 
HS 0837+4717 show no objects within this sphere, corrob- 
orating the classification HS 0837+4717 as a void galaxy. 

4.6. HS 0837+4717 in relation to other XMD BCGs 

One of the goals of detailed studies of the properties of 
individual XMD BCGs is to establish how homogeneous 
is this group of galaxies. Current models suggest various 
evolution scenarios that would result in a large metal de- 
ficiency of a galaxy's ISM. They include both a signifi- 
cant loss of fresh metals via galactic superwinds, related 
to strong starbursts, and/or the infall of unprocessed in- 
tergalactic matter onto normally evolved galaxy. The third 
option is a very slow astration due to very low surface gas 
density and the stabilizing role of the Dark Matter (DM) 
halos (as, e.g., for LSB galaxies). Finally, one of the most 
intriguing options is that of truly young XMD galaxies. 

Known XMD galaxies can, in principle, be the 
products of any of these evolutionary scenarios. 
Therefore, their properties could show significant diver- 
sity. HS 0837+4717 is, in this context, an outlier among 
XMD BCGs and could yield important information on al- 
ternative evolutionary scenarios. 

What are the main differences between HS 0837+4717 
and prototypical XMD BCGs I Zw 18 and SBS 0335-052E, 
both known to show WR features (Sect. I4.3|l . but not an 
excess nitrogen abundance? The BCG discussed here, is 
the most luminous and is probably the largest galaxy in 
this group. The latter implies a rather massive baryon and 
DM configuration, capable in the isolated state of main- 
taining a large fraction of the newly produced metals. Its 
value of M(Hi)/Lb ^0.9 implies either significant pro- 
cessing of neutral gas, or a very strong starburst, or both. 
Also, a significant gas mass fraction could be hidden in 
molecular form. The latter is more characteristic of lumi- 
nous IR galaxies, a class of galaxies tightly connected to 
powerful starbursts due to the merging of gas-rich galax- 
ies. The strong starbursts, in turn, occur when they en- 
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compass a significant part of the galaxy gas mass, what 
is expected in case of strong disturbances. Merging of two 
galaxies provides the strongest possible disturbance. This 
allows the breaking of the gas equilibrium even in very sta- 
ble configurations, and results in the efficient gas sinking 
into the dynamical center of merged object with a subse- 
quent intense SF burst (e.g., Mihos & Hernauist ll996|l . 

The important aspect of star formation in a merger 
is related to the processes taking place in the pre-merger 
phase. As shown by Elmegreen et al. H1993|l . the interstel- 
lar gas of pre-merger galaxy is highly agitated. This results 
in a significant increase of the Jeans mass for forming grav- 
itationally bound molecular clouds (up to 10^-10^ ^q)- 
The latter could be progenitors of the HS 0837-^4717 
very massive starbursts. As discussed, e.g., by Elmegreen 
H1999|l . in regions with very high SF activity one should 
expect a shift of the IMF to higher masses. Similarly, due 
to the lower cooling efficiency at very low metallicity, the 
Jeans mass is expected to be higher, again resulting in 
an upward shift of the respective IMF. This could result 
in a significant increase of the relative fraction of massive 
stars, and a larger effect of WR stars on the local nitro- 
gen enrichment. Thus, the option of merging of two slowly 
evolving XMD LSB galaxies could be a plausible explana- 
tion for the XMD BCGs with such atypical properties. 

4. 7. Summary and Conclusions 

We performed a complex study of one of the very metal- 
deficient BCG HS 0837+4717. The long-slit high S/N ra- 
tio spectrophotometry was complemented by a study of 
the ionized gas velocity distribution, by Hi radio line 21- 
cm total flux measurements and by broad-band imaging 
photometry. The observed properties of this object are 
extreme in several aspects. 

This is the most luminous BCG with a metallicity less 
than or of the order of Zq /20. Its observed unusual colours 
are caused mainly by the very strong [O ill] emission lines 
shifted into the IZ-band due to its significant redshift. 
The factor of 6 nitrogen overabundance in its supergiant 
H II region, in comparison to other very metal-deficient 
BCGs, implies either a non-typical evolution track for HS 
0837+4717, or some very short transitional phase of star 
formation, or both. The ionized gas kinematics appears 
rather disturbed. There is an indication of either counter- 
rotation, or of a supershell. The BCG overall disturbed op- 
tical morphology with the " double- nucleus" central struc- 
ture, and two small appendages at its periphery, may sup- 
port the hypothesis of an advanced merger. Recent images 
of this BCG from the SDSS DRl database corroborate the 
described appearance. 

Among the known BCGs with the well measured N/0, 
there are two objects with the larger N excesses, and three 
more BCGs with an excess of a factor of 3 relative to 
the norm for this type of objects. All of them arc WR 
galaxies; almost all belong to the group of luminous BCGs; 
and all can be interpreted as various stages of mergers. 



This suggests that HS 0837+4717 could presumably be 
the most metal-poor representative of this group. 

From the results and discussion above we draw the 
following conclusions: 

1. The oxygen abundance in HS 0837+4717 is 
12+log(0/H)=7.64. The abundance ratios X/0 
for the elements Ne, Ar, S, CI and Fe are consistent 
with the average values found for low-metallicity 
BCGs by Izotov & Thuan ifT^p . implying the 
primary origin of these elements along with oxygen. 

2. HS 0837+4717 is the most metal-deficient galaxy 
among the luminous BCGs with Mb <-18™0. This 
BCG with its O /H and Mb significantly deviates from 
the Z-L-Q relationship derived for the large BCG sam- 
ple. 

3. We detected the broad emission features characteristic 
of WR stars, and estimated the number of WR stars 
as ~1000. Their ratio to the number of 0-stars is well 
compatible with that predicted by the current models. 

4. Broad low-contrast components of the Ha, H/3, and 
[Qui] AA4959, 5007 emission lines were detected, with 
FWHMs of -1500 km s^^. The fiux ratio of the Ha and 
H/3 broad components implies very strong obscuration 
('^5 mag in B-band) in the emission region or along 
the line of sight. The extinction-corrected luminosity 
of the BCG should be higher than that observed, at 
least, by a factor of 2.5-5. 

5. The position- velocity diagrams for the Ha narrow 
component are not compatible with regular rotation. 
The major axis P-V diagram indicates either counter- 
rotation, or the presence of a giant supershell with a di- 
ameter of —3.5 kpc and an expansion velocity of about 
70 km s^^. Half of the P-V diagram can be interpreted 
as rotation with an amplitude of —50-70 km s^^ at a 
radius of —4 kpc. 

6. We find that this BCG is asymmetric both near the 
center (the inner structure consisting of two compact 
regions —2 kpc apart) and at the periphery. While the 
disturbed morphology of HS 0837+4717 could suggest 
a recent strong interaction, no candidate galaxies are 
found in its vicinity. Along with the very disturbed gas 
kinematics this could be evidence for a hypothesis of 
a recent merging of two gas-rich dwarfs. 

7. The nitrogen-to-oxygen abundance ratio in HS 
0837+4717 is six times higher than that of other XMD 
BCGs. The properties common to HS 0837+4717 and 
five other non-XMD BCGs with the large nitrogen 
abundance excesses, suggest that the large nitrogen 
overabundance could be connected with merger events 
and with a short phase of the related powerful star- 
burst, when many WR stars contribute to the ISM 
enrichment. 

8. HS 0837+4717 is located in a region with very low 
density of galaxies (void), the nearest of them with 
known redshift being situated at —6.4 Mpc h^^. The 
strong isolation of the BCG or its progenitor(s) might 
be the reason for its slow chemical evolution. 
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